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Abstract — Space-Time Block Codes from square complex or- 
thogonal designs (SCOD) have been extensively studied and most 
of the existing SCODs contain large number of zero. The zeros 
in the designs result in high peak-to-average power ratio (PAPR) 
and also impose a severe constraint on hardware implementation 
of the code when turning off some of the transmitting antennas 
whenever a zero is transmitted. Recently, rate | SCODs with 
no zero entry have been reported for 8 transmit antennas. In 
this paper, SCODs with no zero entry for 2 a transmit antennas 
whenever a + 1 is a power of 2, are constructed which includes 
the 8 transmit antennas case as a special case. More generally, 
for arbitrary values of a, explicit construction of 2™ x 2™ 
rate 2±i SCODs with the ratio of number of zero entries to 

the total number of entries equal to 1 2±L2 ll ° 92i ~^ )S is 
reported, whereas for standard known constructions, the ratio is 
1 t±! . The codes presented do not result in increased signaling 
complexity. Simulation results show that the codes constructed in 
this paper outperform the codes using the standard construction 
under peak power constraint while performing the same under 
average power constraint. 

Index Terms — Amicable orthogonal designs, MIMO, orthogo- 
nal designs, PAPR, space-time codes, transmit diversity. 



I. Introduction 

SPACE-TIME Block Codes (STBCs) from Complex Or- 
thogonal Designs (CODs) have been extensively studied 
in Q, 0, 0. 

Let xi,X2,--- ,xt be commuting, real inde terminates. A 
real orthogonal design X of order n and type (<zi, a%, ■ ■ ■ , at), 
denoted as OD(n; a\, a,2, ■ ■ ■ , at) where the coefficients are 
positive integers, is a matrix of order n with entries chosen 
from 0, ±xi, ±X2, • • • , ±cc t , such that X T X = (a-\x\ + 
aix\ + • ■ • + a t Xt )I n where X r denotes the transpose of the 
matrix X and I„ is the n x n identity matrix. 

Amicable orthogonal designs (AODs) are defined using 
two real orthogonal designs of same order but not neces- 
sarily of same type. Let X be an OD(n; m, 112, ■ ■ ■ ,u s ) 
on the real variables x\,X2,--- ,x s and let Y be an 
OD(n; vi,V2, - ■ ■ , ft) on the real variables yx,V2, ■•' , Vt- It is 
said that X and Y are AOD(n; u\, U2, ■■ ■ , u s ;vi,V2, ■ ■ ■ , ft) 
if XY r = YX T . 
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Amicable orthogonal designs have been studied by several 
authors ifTOl . ifTTI to construct complex orthogonal designs. 
The book by Geramita and Seberry fiTI gives a nice intro- 
duction to this topic. 

In the following, we define square complex orthogonal de- 
sign which we use frequently in the rest of the paper. A Square 
Complex Orthogonal Design (SCOD) G(x\,X2, ■■■,Xk) (in 
short G) of size n is an n x n matrix such that: 

• the entries of G(xi,X2, ■■■,£&) are complex linear com- 
binations of the variables x\, Xa, x% and their complex 
conjugates x\ , x% , . .. , x\ , 
. G H G = (|:ci| 2 + ... + |x fe | 2 )I„ where H stands for the 
complex conjugate transpose and I n is the n x n identity 
matrix. 

If the non-zero entries are the indeterminates ±xi,- ■ • ,±Xk 
or their conjugates ±x*, ±x%, ±xJ only (not arbitrary 
complex linear combinations), then G is said to be a restricted 
complex orthogonal design (RCOD). The rate of G is £ 
complex symbols per channel use. 

It is known that the maximum rate 1Z of an n x n RCOD 
is 2±1 where n — 2 a (2b + l),a and b are positive integers 
(21 . Note that the maximal rate does not depend on b. Several 
authors have constructed RCODs for 2 a antennas achieving 
maximal rate 0, (4|, |5|, |6|. In [2|, the following induction 
method is used to construct SCODs for 2 a antennas, a = 
2,3,- • • , starting from 



Gi = 



x l x 2 
X2 Xl 



G„ = 



G a _l — X* a + l l 2 a 

a;„ + il 2 a-i G^i 



(1) 



where G a is a 2° x 2 a complex matrix. Note that G a is a 
RCOD in a+1 complex variables x\, X2, • ■ • , x a +i - Moreover, 
each row and each column of the matrix G a contains only 
a + 1 non-zero elements and all other entries in the same row 
or column are filled with zeros. The fraction of zeros, defined 
as the ratio of the number of zeros to the total number of 
entries in a design, for G a , is 

1 



2 a - a - 1 a 
= 1 - 



= 1-11. 



(2) 



2 a 2 a 
For the constructions in 0, (4), 0, (6) also, the fraction of 
zeros is the same as given by 0. Reducing number of zeros in 
a SCOD for more than 2 transmit antennas (for two antennas, 
the Alamouti code does not have any zeros), is important for 
many reasons, namely improvement in Peak-to- Average Power 
Ratio (PAPR) and also the ease of practical implementation of 
these codes in wireless communication system fl3l . 

For illustration, consider the SCOD G2 of size 4 shown 
below - it is a RCOD, whereas the code Gtjc a ls° shown 
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X 4I + J X 3Q X 3J + J X 4Q x 4/ + J X 3Q x 37 + J X 4Q X 2J + J X 1Q x 2/ + J X 1Q X 1J + 3 X 2Q x 17 + i x 2Q 
x 3J + Jx 4 q x 4J + JX 3 q x 3J + j'x 4 q x 4/ + jx 3 Q x 2J + JXjq x 2J + j'xjq asj^ + jx 2 Q x 17 + jx 2 Q 

x 4 j + J' X 3Q x 37 + J X 4Q x 47 + J X 3Q x 3/ + 3 X 4Q X 1J + 3 X 2Q x l/ + J X 2Q x 27 + J X 1Q x 2/ + 3 X 1Q 
X 3J + J' X 4Q x 47 + J X 3Q X 3J + J X 4Q x 4/ + 3 X 3Q X 1J + 3 X 2Q X 1J + J X 2Q x 27 + J X 1Q x 2/ + 3 X 1Q . 



below, given in (TJ, ifTTj . obtained from Amicable Orthogonal 
Designs, is not a RCOD and there are no zeros in this matrix. 



G, 



Xl 
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Notice that some of the entries of Gtjc can be written as 



(- 



X1—XI+X2—X2) 

2 — 

(lEi— X 2 — X 2 ) _ 

2 ~~ 

„ 2 — 

(Zl+Xj +z 2 — x 2 ) _ 



-(xii - 

-(X2I - 
X2I 

-(xii ^ 



-a; 



™* 

&2, 
-Xl, 



(4) 



i^2Q) = 
•Mq) = 

f JXlQ = 

- j«2Q) = 

where $1 — xu + jx2Q and £2 — X21 + jx\Q are the 
coordinate interleaved variables corresponding to the vari- 
ables xi and X2, where xu and Xiq are the in-phase and 
the quadrature-phase of the variable xi. Single-Symbol ML 
Decodable Designs based on coordinate interleaved variables 
have been studied in |[T2l . For our purposes, it is important 
to note that whenever coordinate interleaving appears, it is 
nothing but a specific complex linear combination of two 
variables, which will have impact in terms of the signaling 
complexity explained subsequently. 
The following code G3 for 8 transmit antennas, 



G, = 



Gy = 
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X2 Xl 


^3 
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X3 xi 
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X2 


-X2 


X3 -X3 




— X4 
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jxi 


3X2 


3X2 


* - * 

3X3 3X3 


3X*4 


3X*4 




-x 2 


X2 


x* 


* 

x* 




-X3 


x 3 


-jx*2~ 


3X1 


3X1 


-3x1 




-3X3- 


-3X3 




-X3 


X3 




x^ 


1 1 


j; 2 


—x 2 


* 

~JX 3 - 


3X*3 




-3x4, 




3X2 


3x2 






X4 


* 
x 3 


—£3 


—^2 X2 


X*l 


xl 


-jXi- 


* 
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3X3 


3x3 
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-3x1 



(5) 



contains 50 per cent of entries zeros. But, Yuen et al, in 
[|7l , have constructed a new rate- 1/2, SCOD ^fc of size 8 
with no zeros in the design matrix using Amicable Complex 
Orthogonal Design (ACOD) ifTTl where Gy is given in (0. 



Observe that for a fixed average power per codeword, due 
to the presence of zeros in G3, the peak power transmission 
in an antenna using G3 will be higher than that of an antenna 
using Gy. Hence, it is clear that the PAPR for the code Gy 
is lower than that of the code G3. Hence, lower the fraction 
of zeros in a code, lower will be the PAPR of the code. In 
[8 1, [9 1, [ 10 1, another rate- 1/2, 8 antenna code with no zero 
entry, denoted by Gtwms shown at the top of this page, has 
been reported. 

Observe that Gtwms has entries that are coordinated in- 
terleaved variables and hence has larger signaling complexity. 

Signaling complexity: 

Notice that some of the entries, for instance xu+jx2Q and 
X21+3X1Q, in Gtjc and Gtwms are co-ordinate interleaved 
versions of the variables xi and X2- Suppose the variables 
xi and X2 take values from a regular (rectangular) 16-QAM 
rotated by an angle 9. Though rotation does not affect the 
full-diversity of the code, the coding gain depends on 9 and 
hence non-zero value of 9 may be desired. Now the antenna 
transmitting xi chooses one of the 16 complex numbers for 
transmission whereas the antenna transmitting x\j + jx 2 Q 
will be choosing one of 16 x 16 complex numbers since 
the components xu and X2Q take independently 16 values 
each. This will increase the number of quantization levels 
needed in a digital implementation for signals transmitted 
in this antenna. We will henceforth refer to the number of 
quantization levels needed in such a digital implementation as 
"signaling complexity". Notice that designs which have entries 
that are linear combinations of several variables increase the 
signaling complexity of the design. Accordingly, the signaling 
complexity of G2 given in (0 is less than that of the code on 
the right hand side of (0. Similarly, the signaling complexity 
of Gtwms is larger than that of Gy. 

Notice that by multiplying the matrix G3 with a unitary 
matrix, the resulting matrix will continue to be a SCOD with 
different number of zeros and it is not difficult to find unitary 
matrices that will result in a design with no zero entries. 
However, such a design is likely to have a large signaling 
complexity which needs to be avoided. Obtaining a unitary 
matrix which reduces the number of zero entries while not 
increasing the signaling complexity is a nontrivial task which 
is the subject matter of this paper. 

In this paper, we provide a general procedure to construct 
SCODs with fewer number of zeros compared to known 
constructions for any power of two number of antennas 
(greater than 4), without increasing the signaling complexity. 
Our contributions are summarized as follows: 

• Maximal-rate SCODs with no zero entry and minimum 
signaling complexity for 2 a transmit antennas whenever 
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a + 1 is a power of 2, are constructed which includes the 
8 transmit antennas case as a special case. This matches 
with the construction given in Q for 8 transmit antennas 
and beats the codes in (8), (9), iTTOl for 8 transmit 
antennas in terms of signaling complexity. 

• More generally, for arbitrary values of a, explicit con- 
struction of 2 a x 2 Q , rate ^±1 SCODs with the ratio of 
number of zero entries to the total number of entries equal 
to 1 - 2±1 2 l/ ° 92( ^ft ) J is reported. Note that when a + 1 
is a power of two, our codes have no zeros. When a + 1 
is not a power of two, it is conjectured that SCODs with 
smaller fraction of zero entries with rate and same 
signaling complexity do not exist. Our construction gives 
fewer number of zero entries compared to the well known 
constructions in 0, H, 0, 0. 

« Our construction is based on simple premultiplication of 
the code in (HJ by a scaled unitary matrix consisting 
of only +1,-1 or 0, whereas the constructions in Q, 

10, HD, ifTUl depend on the existence and availability of 

AODs urn . 

• A general procedure to obtain the scaled unitary matrix 
that leads to a SCOD with small number of zero entries 
is given. 

• It is shown that the new codes presented in this paper 
admit a recursive relation similar to that admitted by G a . 

The remaining content of the paper is organized as follows: 
In Section [II] we prove the main result of the paper given by 
Theorem Q] In Section |IIIJ we give a procedure to compute 
the premultiplying matrix using which we can get the SCODs 
of this paper straightaway from the well-known construction 
given by ([T). The PAPR of the new codes constructed is 
discussed in Section[lV] Simulation results are given in Section 
[Vl A brief summary and a conjecture constitute Section [VTJ 

11. Construction of SCODs with Low PAPR 

SCODs given in |2] contain a large number of zeros and 
the fraction of zeros in the code increases as the number of 
transmit antenna increases. Note that these codes are RCODs 
and hence of least decoding complexity as well as least 
signaling complexity. It is possible to obtain an orthogonal 
matrix with fewer zero, if we premultiply and/or post-multiply 
the given orthogonal design matrix by some unitary matrix, but 
the resulting orthogonal design need not be a RCOD. So care 
must be taken in how we choose these premultiplying or post- 
multiplying matrices such that the code obtained after applying 
these matrices, does not contain complex linear combination 
of the variables which will increase the signaling complexity. 

There exists a unitary matrix which when pre-multiplies the 
code G3 obtain a code which contains no zero in the matrix 
and none of the entries of this new code is a complex linear 
combination of variables and thus the signaling complexity 
is not increased. The unitary matrix corresponding to G3 is 
-i=Q( 3 ) where Q^ 3 ) is given by the matrix on the left hand 
side of Here —1 is represented by simply the minus sign 
(throughout the paper) and the resulting no zero entry SCOD 
is H 3 where the matrix V2H3 is shown on the right hand side 
of ©. 
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1000 1 
0100 1 
0010 1 
0001 1 
0001- 
0010 0- 
0100 0-0 
_ 1000 0- 

Towards identifying such premultiplying matrices for the 
general case, we label the rows of G a as Rq, Ri, ■ ■ ■ , R 2 a-i- 
The column index also varies from to 2 a - 1. Let be the 
set of column indices of the non-zero entries of the i-th row 
Ri of the matrix G a . The following lemma describes 
for alH = to 2 a - 1. 

Lemma 1: Let a be a positive integer and G a be a COD of 
size 2 a x 2 a in (a+1) complex variables x\, ■ ■ ■ , x a +i as given 
in (HJ. Let i be a positive integer between and 2 a — 1. Let 
the radix-2 representation of i be (i Q -i,i a -2, • •■ , *o) where 
i a -i is the most significant bit. Then 



(6) 



N 



(a) 



{»} U {i + (-1)*^ I j 



0,' 



1} 



or equivalently, N-"' = {i}U{i®2 J \ j = to a-1} where © 
denotes the component-wise module 2 addition of the radix-2 
representation vectors. 

Proof: The proof is by induction on a. The case a = 1, 
corresponds to the Alamouti code Gi. We note that — 
{0, 1} and N± = {0, 1} as given by the expression of 
for i = 0, 1. So for a = 1, the lemma is true. Let the lemma 
be true for all a < n. Then, we have 



= {i}U{i+(-l)*'2i\j = 0, 



1} (7) 



for all i = 0, 1, • • • , 2" — 1 and we need to prove that 



N, 



(n+l) 



= {i}U + 2*'| j = 0,- 



(8) 



for all i = 0, 1, • ■ ■ , 2 n+1 — 1. For a = n + 1, we have the 
radix-2 representation, i = (i n ,i n _i, ■ ■ ■ ,i Q ) and 



G n 

h2l2™ 



b n+2 



I 2 n 



G 



H 



(9) 



We have the following two cases: 

Case (i) < i < 2" — 1 : In this case i n = and the term 
i + (— l) ln 2 n in ^ corresponds to the non-zero location in 
the —x* n+2 \ 2 n part of G„ + i and the nonzero locations in the 
G„ part is given by the remaining elements of ([8]) which is 
nothing but N^ n) . 

Case (ii) 2" < % < 2 n+1 - 1 : In this case i n = 1 for all 
values of i in the range under consideration. Then, the term 
corresponding to j = n in ([8]l is z — 2™ which corresponds to 
the non-zero term in the x n + 2 I 2 n part of the matrix ([9j. Also, 
every term of the form i + (— l) Jj 2 J ; j — 0, 1, • • • , n— 1 in ([8]) 
will be same as a term in (0 with 2™ added to it. This takes 
into account all the non-zero entries in the G^ part of (0. ■ 

Example 2.1: In this example we compute the sets for 
a = 2 and 3. For a = 2, the possible values of i are 0, 1, 2 
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TABLE I 
M a AND M' a FOR a = 3, • • ■ 9 



a 


3 


4 


5 


6 


7 


8 


9 


M a 


{3} 


{3} 


{3,5} 


{3,5,6} 


{3,5,6,7} 


{3,5,6,7} 


{3,5,6,7,9} 


M'a 


{7} 


{7} 


{7,25} 


{7,25,42} 


{7,25,42,75} 


{7,25,42,75} 


{7,25,42,75,385} 


d 


2 


3 


3 


3 


3 


4 


4 



and 3, while for a — 3, i takes value between and 7. 

A^ 2) = {0} U {0 ® 2°, ffi 2 1 } = {0, 1, 2}, 



(«) 



N 
N. 
N. 



(2) 
1 

(2) 
2 

(2) 



{1}U{1©2°,1©2 1 } = {1,0,3}, 
{2} U {2 ©2°, 2© 2 X } = {2, 3,0}, 
= {3} U {3 © 2°, 3 © 2 1 } = {3, 2, 1}. 

={0}U{1,2,4} = {0,1,2,4}, 
N[ 3) = {1}U{0,3,5} = {1,0,3,5}, 
= {2} U {3, 0,6} = {2, 3, 0,6}, 



AT, 

A^ 



(3) 
1 

(3) 
(3) 

r(3) 



{3} U {2, 1,7} = {3, 2, 1,7}, 
{4} U {5, 6,0} = {4, 5, 6,0}, 
{5} U {4, 7,1} = {5, 4, 7,1}, 
{6} U {7, 4, 2} = {6, 7, 4, 2}, 



Similarly, x = joix = j® 2 k ' 2 for some < fc 2 < a — 1, as 

But if x = i, then x ^ j, as i ^ j. So, we have following 
three cases: 

(i) a; = i and x = j © 2 fc2 , 

(ii) x = i © 2 fcl and x = j, 

(iii) a; = i © 2 fel and x = j © 2 fe2 , fc x ^ fc 2 (as i 7^ j). 

For the case (i) & (ii), we have i = j © 2 k2 & i © 2 fel = j 
respectively and in both cases, dist(i,j) = 1. For the case (iii), 
we have i © 2 fcl = j © 2 fea , which means the dist(i,j) = 2. 
So MHD of S is less than 3. ■ 
For a given a, let d be the positive integer such that 2 d ~ 1 < 

a < 2 d and a = Yl'j-a a fi j •> a o e F 2- Note tnat a d-i = 1- 
Define 



M a = {0 < x < a \ x jt 2 k for any fc = 0, 1, • ■ • } (10) 



and 



M. 



; = /2 x - 1 + ^a: i 2 2i - 1 | x = J2 x i 2 ' &M a ,Xj eF 2 }. (11) 



3=0 



j=0 



A r r={7}U{6,5 ! 3} = {7,6 ! 5,3}. 

(3) (3) 

Notice that N: n Af ■ = if i © j = 7, where represents 
the empty set. 

Definition 1: Two rows Ri,Rj of G a are said to be non- 
intersecting if N^ a) n A r j Q) = 0. 

The following lemma is needed to prove Lemma|4]which in 
turn is used in the proof of the main result given in Theorem 

m 

Lemma 2: Let V a denote the set of all radix-2 representa- 
tion vectors of the elements of the set {0, 1 , • • • , 2 a — 1} and S 
be a subset of V a . Then N^ a) n A/ J (a) = for all i,jeS,i^= j 
if and only if the minimum Hamming distance (MHD) of S 
is greater than or equal to 3. 

Proof: We first show that A^ a) n A"j a) = <\> for all i,j £ 
S,i 7^ j implies that the MHD of S is greater than or equal to 
3. Equivalently, if the MHD of 5 is 1 or 2, then there exists 
i,j e S,i ^ j, such that A^ a) n N^ a) ^ (/>. Assume that the 
MHD of S is 1 or 2, then their exists i,j G S, i j, such 
that dist(i,j) = 1 or 2. So, either % = j © 2 k for some 
k e {0, 1, ■ • • ,a - 1} if dist(i,j) = 1, or i = j ® 2 kl © 2 fe2 
for some ki,k2 € {0, 1, • • • , a — 1}, fci 7^ fc 2 if dist(i,j) = 2. 
In the first case, i S n AT, ■ a ' and in the second case, 

(i © 2 fel ) e A^ a) n as i © 2 fel = j © 2 fe2 . 
For both cases, A^ (a) n Afj a) ^ 0. 
Next we prove that if the MHD of S is at least 3, then N^ a) n a ~ d - Notice that 2*- 1 7^ 2 2 "- 1 for j = 0, 1, • • • , d - 1 as 

M Q) = for all i,jeS;i^ j, or equivalently, if for some x * V ■ So ' wt{ * x ') = 1 + wt W for 4111 s e M « where wt ( x ) 

J , f„i f„i , stands for the Hamming weight of x. Now wt(x) > 2 as x is 

i,jeS,i^ j, m a) n A^] Q) 7^ 0, then MHD of 5 is less than . - „ „ B . f n * „ V y " 

' ' ■ not a power of 2. So wt{x j > 3. 

Similarly, wt(x'®y') = 2+?«t(xffiy) for all x, y € M OJ x 7^ 

y. Now wt(x © y) > 1 as x 7^ y, which implies that wt(x' © 

y') > 3 for all x', y' € M' a . 

In general, wi^ffi^ffi- • -ffiy^.) = k+wt(yi®y 2 ®- ■ -®yk) 

for k < a — d, y[ 7^ y' 2 7^ • • • 7^ y^.. So for all fe > 3 and 

k < a — d, wt(y[ ® y' 2 ® ■ ■ ■ ® y' k ) > 3. Now there exists 



Note that the number of elements in M a is a — d. Moreover, 
M a C Mb and M' a C whenever a < b. When a is a power 
of 2, M a = M a -i and - M^. 

Example 2.2: The sets M a and Af^ for a — 3 to 9 are 
shown in Table U at the top of this page. 

Lemma 3: Let M' a be as defined in (TTTb . View as a 
subset of V = ¥2 by identifying each element of M' a with its 
radix-2 representation vector of length a. Then the MHD of 
the linear space spanned by M' a , denoted by S, is 3. 
Proof: The subspace S C V is given by 

ct rv^a- d— 1 / 

5 = {2^=0 c iVj 
0, 1, • • • , a — d — 1. 

Observe that the map given by 



j/j e M' a } where Cj e F 2 for j 



f : M a -» 

= E?i 1 ^^ a / = 2*- 1 + E J to 1 



2 J - 



(12) 



is one-one. Thus, the size of M^, denoted as \M a \ is also 



3. 

Let i,j £ S and i 7^ j. We have 



(o) 

r(o) 



{i © 2 fe |fc = 0, • • ■ , a - 1} U {i} and 



JV) o; = {i €D 2 & |/c = 0, • • - ,a-l}U{j}. 



As W (a) n W (a) 7^ 0, let x € A^ a) n AT (a) . 

We have j; = i or x = i 2 1 for some < fci < a — 1, as 
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an element in S, for instance 7, whose Hamming weight is 3. 
Thus, the MHD of S is 3. ■ 

Lemma 4: Let a and d be non-zero positive integers such 
that 2 d ~ 1 <a<2 d and V a = {0, 1, • • • 2 a - 1}. Then, there 
exists a partition of V a into 2 d subsets c\ a ^ , j = 0, 1, • • • , 2 d — 
1 each containing 2 a elements, such that for any two distinct 



elements x, y G cj a) , j G {0, 1, • • • , 2 d - 1}, we have N^ a> n 

Proof: We identify the set V a with Ff by viewing each 
element of V a with its radix-2 representation vector. Let M' a be 
as given by (fTTT i and S be the sub-space of V a spanned by the 
radix-2 representation vectors (of length a) of the elements 
of the set M' a . The number of elements in S is 2 a ~ d . By 
Lemma [3] the MHD of 5* is 3. Now we define a relation ' ~' 
on V a as follows: For all a, b € V a , a ~ b, if a © b G 5. We 
observe that this is an equivalence relation as for all a, 6 and 

C G V a , 

1) a~aasa©a = 0eS, 

2) a ~ i =J- 6 ~ a as a ® & e S, implies that & © a G 5. 

3) a ~ & and 6 ~ c, then a ^ c, as a ffi 5 e S and 6 © c G 
S 1 , together imply a © c G 5. 

The number of equivalence classes is t^tj- — 2 d and these 



(a), 



1. 



equivalence classes are denoted as C\ , z = 0, 1, • ■ • , 2 d - 
For any one equivalence class , the elements in cf"' 1 are 
given by {x © s|s € 5} for some a; G Cf . Now the MHD 
of the class c\ a \ is also equal to the MHD of S which is 3. 
By lemma 12] N ( x a) n N^ a) = for all x, y G C*| a) , i = to 
2 d -l. ■ 
The following example illustrates the partition of V a into 
the subsets c\ a \ i = to 2 d - 1, for a = 3, 4, 5 and 6. 

Example 2.3: (i) Let a = 3. V3 is partitioned into 4 
classes Cp 3 ', Cj; 3 '', and Cg 3 ^, each containing 2 elements. 
We have already seen that M3 = {3} and M3 = {7}. 
Cf a) = {i, i © 7} for i = 0, 1, 2 and 3. Explicitly, 



Ci 



(3) 



{0, 7}, C7p } = {1, 6}, = {2, 5}, = {3, 4}. 



-,(3) 



(3) 



(ii) For a = 4, 

C7( 4) ={0,7},C7{ 4) = {l,6},C7( 4) -{2,5}, 



:*) 
-;<4) 



{3, 4}, a 



•4) 



={10,13} ) C( 4) = {11,12} 
(iii) For a = 5, 



{8,15},C7f 



{9,14}, 



C^ 5) = {0, 7, 25, 30}, CT = {1,6, 24, 31}, 



(5) 



a 



(5) _ 



{2,5,27,28},C| 5) _={3,4,26,29}, 



CY ] = {8, 15, 17, 22}, C| 5) = {9, 14, 16, 23}, 



a 



{10, 13, 19, 20}, Cf ] = {11, 12, 18, 21}. 



(iv) For a = 6, 
C<, 6) = {0, 7, 25, 30, 42, 45, 51, 52}, 



(6) 

h 



{1,6,24,31,43,44,50,53}, 
{2,5,27,28,40,47,49,54}, 



C;T = {3, 4, 26, 29, 41, 46, 48, 55}, 
c] 6) = {8, 15, 17, 22, 34, 37, 59, 60}, 



(6) 

h 
.(«) 



{9,14,16,23,35,36,58,61}, 
{10,13,19,20,32,39,57, 62}, 



C [ 7 °> = {11, 12, 18, 21, 33, 38, 56, 63}. 

Theorem 1: Let a and d be non-zero positive integers 



and 2 d - x < a < 2 d . There exists a SCOD H a of 
size 2 a x 2 a with entries of the matrix 2t"H„ consisting 
±X\,±X2, • • • , ±x Q +i or their conjugates, such that the code 
has rate 1Z = and the ratio of number of zeros to the total 
number of entries of the matrix is equal to 1 — 72. ■ 2 L lo S2 tz J . 

Proof: The SCOD H a satisfying the required rate and 
fraction of zeros in the matrix is obtained from the COD G a 
of size 2 a x 2 a (given in ((TJ) as follows: 

The rate TZ of the COD G Q is 3±1. Using Lemma 14] 2 a 
rows of the COD G a can be partitioned into 2 d groups with 
each group containing 2 a ~~ d rows such that any two distinct 
rows from any of 2 d groups, is non-intersecting. If we add or 
subtract all the rows in a given class, the resulting row will not 
have any entry which is a linear combination of two or more 



variables. Labeling the groups as C^' , c[ a ^ , 



,G 



(a) 



we 



define the 2 a d x 2 a matrices formed by the rows of G a 
which are in for i = to 2 d — 1. Now form the matrix 



B' 



B„ 
Bi 



(13) 



The matrix B' is of size 2 a x 2 a and it is related to G a by B' = 
PG a where P is a permutation matrix of size 2 a x 2 a . We 
consider a Hadamard matrix H of size 2 a ~ d x 2 a ~ d containing 
1 and -1 such that H r H = 2 a -%a- d . Let B 4 = HB,. The 
required matrix H a is 



B 
Bi 



Bo 



(14) 



For two matrices A = [a%j] and B, the tensor product of A 
with B, denoted by A ® B, is the matrix [ojjB], Let H = 
I 2 d O H. We can write H a = 2-^HB' = 2-^HPG a . 
Now H Q is a SCOD if and only if 2 ~ HP is an unitary 
matrix. As P is permutation matrix, it is enough to prove 
that 2 2~H is an unitary matrix. Indeed, H r H = I 2 d <E> 
(H r H) = 2 Q - d I 2 a, thus H a is a SCOD. 

The number of locations containing in any row of H Q is 
2 a — (a + 1)2°^. Hence the fraction of zeros in H a is equal 

2°-(a+l)2°-' i _ j_ a±L 2 a - d . NOW 2 a - d < < 2 a - d+1 

2 a — a+1 



to 



as 2 d > a + 1 > 2 d ~ 1 . So, a - d < logo 4r < a - d + 1 and 
a-d = [log 2 ^jj = [log 2 ij. 

Thus the fraction of zeros is 1 - TZ ■ 2^- loS2 tcJ . ■ 
The proof of Theorem [T] suggests a recipe to construct 

the SCOD H a with the fraction of zeros specified in the 

statement, from a COD G a given in tjl). The following 

example illustrates this recipe. 

Example 2.4: We consider the construction of rate-1/2, 

8x8 COD with no zero in the matrix. Following the recipe 

described above, the permutation matrix P and H are given 

by 
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- 10000000- 




- 1 


10 


00 


00 


o-i 


00000001 




1- 


-10 


00 


00 





01000000 







01 


10 


00 





00000010 


,H = 





01- 


-10 


00 





00100000 





00 


01 


10 





00000100 







00 


01- 


-10 





00010000 







00 


00 


01 


1 


. 00001000_ 




. 


00 


01 


01- 


-1. 



respectively. 


The matrix H3 




2 2HPG3 is 


given 


by 




XI 


-x\ 


-X*3 


X4 


— x% 


—xa 


X2 


xl 




Xl 




-x*3 - 


X4 


—xl 


X3 


—X2 


-xl 




X2 


X* 


X4 — 


X*3 


-xa 


—xl 


Xl 


~X*2 


1 


X2 


X* 


— X4 — 


X*3 


X3 


—xl 


-Xl 


X*2 




X3 


X4, 




X*2 


-X2 


Xl 


—xl 


X*3 




X3 


—X4 


xl 


X*2 


X 2 


-Xl 


X4 


~X*3 




Xi 


X3 


-X 2 


Xl 


Xl 


X*2 


X3 


—xl 




— Xi 


X3 


-X2 


Xl 


-xl 


~X*2 


-X3 


—xl 



which is a row permuted version of the code given in ©. 

III. Premultiplying matrix 

In this section, we present a procedure to compute a matrix 
denoted as of size 2 a x 2 a which when pre-multiplies G a , 
along with an appropriate scaling factor, the resulting matrix 
H a is the SCOD with desired fraction of zeros. The scaling 
factor, when multiplied to the matrix 0/ a ', makes it a unitary 
matrix. 

In order to construct the matrix Q( a ', we first associate a 
2 X x 2 X matrix Q x to each x in M a . The (i,j')th element of 
Q x , denoted by Q x (i,j), is defined as follows: 

For i = to 2 X - 1 - 1, 

1) Q x (m) = 1; 

2) Q x (i,i®x') = 1; 

3) Q x (i®x',i) = 1; 

4) Qx(i®x',i®x') = -1; 

5) Qx(i,j) = for all other values of i and j; 

where x' is given by ( fT2b . Define a 2 a x 2 a matrix Q x as 
Q x — I^a-x Q x where I^a-x is the identity matrix of size 

i^a—x ^a—x 

Example 3.1: In this example, we compute Q x for x £ M4. 
We have M4 = {3}. The matrix Q3 is the matrix shown on 
the left hand side of (0 without the scaling factor 

Towards the construction of premultiplying matrix, we need 
the following result which says that Q x and Q y commute for 
all x,y&M a . 

Lemma 5: Q x Q y — Qj,Qa; for all x,y G Ma- 
Proof: Let the ith row of Q x be Q* . There are 2 a rows 
and 2 a columns for the matrix Q x with non-zero entries either 
1 or -1. 

Fix i. Let Q x (i, I) = a, q e {0, 1,-1} for I = to 2 a - 1. 
So Qj. = (c , ci, • • • , c 2 a_i). We write (c , c 1( • • • , c 2 a^ 1 ) as 
J2i—o 1 Q2 ! and this correspondence is unique in the sense that 
no two distinct vectors will produce the same value under the 
above correspondence. Then, we have Q^, = Y2i=o 1 c '^- Let 
the radix-2 representation of i and j be (i a -i,i a -2, ••■ , io) 
and (ja-iJa-2, ■■■ , Jo) respectively. 

Note that = (-l)*—^ + 2 1 ® X ' for i = to 2 a - 1. 



Moreover, Q x is a symmetric matrix for all x <E M a . Let 
K = Qn-Qj, and the (i,j)th entry of K be K(i, j). It follows 
that 



K(i,j) = <^ 



U 




It I t 


B J 


/T* ffi r*J nJ nr-t 1 

¥ {U,x ,y ,x f 


By) 


{- 


1 \i-r-l 
■1J * 1 


It I k 


B J 


= V 




1 

\~ 




It I ( 


B J 


— X 




1 
1 




if i t 


v J 


— r r f rTi 1/ 




,«- 


^iyx-l+jy-1 


if * € 


BJ 


= , 









if * € 


BJ 


£ {0,x',y',x'( 


By'} 


(- 




if z € 


B j 


— y' 




(- 


-l)*y-i 


if ii 


B j 


= x> 




1 




if z € 


B j 


= x'®y' 






.\\3x-l+iy-l 


if i € 


BJ 


= . 





K(j,i) 



Let a; = X^f^o 1 ^ 2 ' and 2/ = Sf^^ 2 '- We have y' = 
2J/-1 + j/j2 2l_1 . Note that the (x - l)th component of 

y', i.e., the coefficient of 2 X ~ 1 in the radix-2 representation 
of y' is zero as x ^ y and a; is not a power of 2. Thus 
i x -i = jx-i when i G) j = y' ■ Similarly j y -\ — i y -i when 
i®j — x 1 and i x -\+j y -\ = j x -\+i y -i when i©j = 0, i.e., 
i = j. Thus K(i, j) = K(j, i) for all i, j e {0,_1, - ■ , 2 a - 1}. 
So, K = Q K Q y is symmetric. Then, K r = Q^Q^ = K = 
QxQy As Q y and Q x are symmetric matrix, they commute. 

■ 

Let Q( a ' = \[ Q x , which is well defined since the 

xeM a 

product of these matrices does not depend on the order these 
matrices are multiplied. The following theorem asserts that 
Q( Q ' so constructed, will produce a SCOD with the desired 
fraction of zeros. 

Theorem 2: Let a and d be non-zero positive integers and 
2 d ~ 1 < a < 2 d . Then H a = 2-^Q( a )G a is the desired 
SCOD with the rate and the fraction of zeros as specified in 
Theorem [T] 

Proof: We have to prove that H a = 2~ V^Q^Ga is a 
SCOD of size 2 a x 2 a with rate K = &± and the fraction 
of zeros is equal to 1 — 72. - 2L'°S2 ttJ. Since we obtain H 
by pre-multiplying a constant matrix to G a , the rate remains 
same. Thus, it is enough to prove that H a is a SCOD and it 
contains the desired fraction of zeros. As G a is a COD, so 
H a is a COD if 2 j-QW is an unitary matrix. Moreover, if 
each row of H a contains 2 a ~ d (a + 1) non-zero entries, then it 
contains the required fraction of zeros. It is easy to note that 
if the matrix Q^ a ^ has 2 a ~ d non-zero elements in each of its 
rows, then each row of H a will have 2 a ~ d (a + 1) non-zero 
entries. The column co-ordinates of the non-zero entries of zth 
row of Q( a ) be such that the resulting matrix H a will not have 
any entry which is linear combination of complex variables, 
i.e., those rows will be added or subtracted which possess 
non-intersecting property. Thus we have to prove following 
two claims: 

1) 2 i _ Q(°) is an unitary matrix; 

2) The column co-ordinates of non-zero entries on the i- 
th row of O/"' is given by the set Sf = {s © i\s S 
S a } where the subspace S a C Ff is spanned by the set 
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We know that \M a '\ is a-d, thus \S a \ is 2 a - d . 
First we prove claim 1). Let x £ M a . Consider the i-th 
and j-th column of Q x which are denoted as Q l x and Q£ 
respectively. The inner product of Q' x and Q^, denoted as 
(QxjQi) is as follows: 

1) for j^i,i®x', (C&Qj) =0, 

2) for j = i®x', (Q*,Qj) = 1-1 = 0, 

3) FQri = t,<Q*,Q£) = l + l = 2. 

So Q = 2I 2 x where I 2 * is an 2^ x 2 a: identity matrix. 
Now Q^Q r = I 2 a-x © Q^"Q X which implies that Q X Q X = 
2I 2 a and Q( a ) T Q( Q ) = 2 a ~ d I 2 a. Thus, 2"^Q( a ) is an 
unitary matrix. 

Next, we jirove claim 2) by induction on a. Let a = 3, then 
Q( 3 ) = Q 3 = Q 3 . The i-th row of Q3 contains non-zero 
entries only at i and % © 7, corresponding to the elements of 
{s © i|s G S 3 } where S* 3 = {0,7}. Observe that S 3 is the 
subspace spanned by M' 5 — {7} C F 2 . 

Let it be true for a < (n — 1). Under this assumption, the col- 
umn co-ordinate of the non-zero entries on the ith row of the 
matrix Q( n_1 ) is given by the set = {s © i\s € S™ -1 } 

for i = to 2™ _1 — 1. We have the following two cases: 
Case (i) n is a power of 2: In this case M' n = M' n _ x 
and = (I 2 © Q*"- 1 )). fact 2) is trivially satisfied. 

Case (ii) n is not a power of 2: M n = M„_i U {n} and 
Q(") = Q„(I 2 © QC"- 1 )) = Q„(I 2 © Qt™" 1 )). The i-th 
row of Q„ is given by = (- l)*"- 1 2* + 2 2 ®"' for i = 
to 2" - 1 and j-th row of Q^"" 1 ), denoted by Q("- 1 )J is 
X] s gs" _1 (~ 1)" S 2 S where a s is either +1 or —1, depending 
on s. The (i,j)-th entry of QW is given by the inner product 
of the i-th row of Q„ with j th row of (I 2 © Q^"" 1 ') as 
q("-i) an( j (I 2 ©Q( Tl-1 )) both are symmetric matrix. Viewing 
S n-i c ¥ n-i^ as a subset of by identifying F™" 1 as a 
subspace of F 2 , it is possible to express S" in term of S™^ 1 . 
The column co-ordinate of the non-zero entries at the i-th row 
of Q(») is S? = 5J'- 1 U{«ffin / |a G for i = to 2 n -l. 

Thus = {s®i\s G 5 n }. ■ 

The following theorem shows that the matrix H a with re- 
quired rate and the fraction of zeros as specified in TheoremQ] 
can also be constructed recursively, in a similar fashion as for 
G a , and it is done using the premultiplying matrix QW. 

Theorem 3: Let a be a non-zero positive integer and H a be 
the SCOD as stated in Theorem |2] Then H a+ i is constructed 
recursively using H a as follows: 



H 



0+1 



H' a+1 when (a + 1) is a power of 2 

^Qn+iH'a+i otherwise , 



where 



H' 



0+1 



H a (xi,X2, ■ ■ ■ ,x a+ t) 
z Q+2 H a (l,0,--- ,0) 

Proof: We have 



—x 



a+2 ±1 -a\^, 0. 



-x 2 , ■ ■ ■ 



•••,o) 

-Xa+l) 



M a+ i = 



M a 



when (a + 1) is a power of 2 , 



M a Li{a+l} otherwise 



Now Q( Q) = J] where Qz is 2° x 2a matrix. 

xeMa 

To indicate the dependence of the size of Q x on a, we write 



We have Q a x = I 2 a-«> © and QJJ+ 1 = I 2 a+i- 
QS +1 = I 2 © Q%. 

Moreover, Q( a ) = J] QS an d Q (a+1) = II Q* +1 - 

16M, xeM a+ i 
Observe that, if a + 1 is not a power of 2, then 



Q x . So, 



Q 



(o+l) 



Q^-(l 2 ©Q (a) ). 



Since Q"tJ — Qa+i> we nave 



Q 



Q 



(„.) 



(a+l) = J A 2 

1 Qa+i • (I2 © Q (a) ) otherwise 



when (a + l) is a power of 2 



From Theorem [2] we have H a = 2 2 Q( a )G a where d 
is given by 2 d ~ 1 < a <2 d . Hence, 



H a+ i = 
Let 

Then, 
H a +i 



2 ~s~ ■ (I 3 ® Q [a> ) ■ Ga+i when (a + 1) is a power of 2, 
2 _o-d+i _ ^^^^ _ ^ ^ ^ {a) ^ _ G ^^^ otherwise. 

H; +1 = 2-^(I 2 ©Q( a ))G a+1 . (15) 

H^ +1 when (a + 1) is a power of 2 , 

^=Qa+iH' a+1 otherwise. 



H' a , j is constructed using Ff a as follows: 
We have from (Q}, 



G a +i 



(16) 



G a -< +2 I 2 <* 
^0+212° G^- 

From the construction of G a and H a , it follows that 

I a . = G.(l, 0, ■ • ■ , 0), Q (a) = 2^H a (l, 0, ■ • • , 0), 
G^ = G a (xi,X2, ■ ■ ■ ,x a +i) = G (ar*, -X2, ■ • • , —x a +i). (17) 
Using ( fT5l ), ( [16] ) and ( fT7| ), we have 



For a 
follows: 



Hl + l(ll,l2, • • • ,x a+2 ) 

H a (xi,x-»,'-' ,x a+1 ) -x* +2 H a (l, 0, • • • ,0) 

X a+2 H (1, 0, • • • ,0) H- a (x*, — X 2 , ■ ■ ■ ,-Xa+l) 

3, Q^ 3 ' is given in ([6j. For a = 4, O/ 4 ) is as 



1000 
0100 
0010 
0001 
0001 
0010 
0100 
1000 
0000 
0000 
0000 
0000 
0000 
0000 
0000 
0000 




1 
10 
10 
-00 
0- 
0- 












10000 
00000 
00000 
00000 
00000 
00000 
00000 
-0000 
01000 
00100 
00010 
00001 
00001 
00010 
00100 
01000 



oooo- 








1 
10 
10 
10 
-000 
0-0 
0-0 
0-. 



The SCOD obtained by premultiplying G 4 with 2~1?Q (4) , 
is shown on the left hand side at the top of the next page. 
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TABLE II 

COMPARISON OF POWER DISTRIBUTION CHARACTERISTICS 





16 Tx; QPSK 


16 Tx; 16 QAM 


32 Tx; QPSK 


32 Tx; 16 QAM 




Peak/ave 


Po 


Peak/ave 


Po 


Peak/ave 


Po 


Peak/ave 


P) 


SCODs G a 


3.2 


0.6875 


11.52 


0.6875 


5.33 


0.8125 


19.2 


0.8125 


Codes in this paper 


1.6 


0.375 


5.76 


0.375 


1.33 


0.25 


4.8 


0.25 



TABLE III 

Variation of fraction of zeros with the number of antennas 



a 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 


15 


16 


/«(H«) 





3 


•2 


i 







B 




4 


3 











/.(Ga) 


1/2 


-Sl- 
id 


1fi 


-§7- 

fi4 


L'2l) 
12R 


2Kfi 


R12 


1 024 


204S 


40flfi 


-8T7S- 

R1S2 


IOmbh 

Ifl IS 1 


327fiR 


§§5iy 

RfiKSfi 



V2 



For comparison, on the right hand side, we have displayed the 
SCOD G4 to compare the number of zeros. The premultiply- 
ing matrix Q( 5 ) for 32 antennas corresponding to a — 5 is 
displayed in Fig. [5] and the resulting code H 5 in Fig. [6] 

IV. PAPR OF THE NEW CODES 

In Table HU the SCODs constructed in this paper are 
compared to the known SCODs given in (flj. For some fixed 
number of transmit antennas, we observe that the PAPR of 
these new codes is less compared to the known SCODs. 
In fact, it is easily seen that as the number of transmit 
antenna increases, these codes outperform the existing CODs 
significantly as far as PAPR is concerned. Quantitatively, for 
QAM signal set, the PAPR of a SCOD for 2 a antennas given 
by CD, is -ttt> while it is ^ — w — for the SCOD of 

same size constructed in this paper. The codes constructed in 
this paper contain fewer zeros than the well-known SCODs. 
Hence, the probability Po that an antenna transmits a zero 
symbol (or switched off) is less in these codes compared to 
the codes given in ((T). 

Another interesting fact to note is that the new SCODs for 
2 a transmit antennas, contains no zero entry when a + 1 is a 
power of 2. For all other values of a, the fraction of zeros keeps 
reducing starting from a = 2 l to 2 l+1 — 1, I a positive integer. 
On the other hand, the fraction of zeros keeps increasing as 
a (> 3) increases for the codes given in (|T). Table [Til] shows 
the variation in fraction of zeros (denoted as f z ) for proposed 
codes H a and the SCODs G a for a = 3 to 16. 

V. Simulation Results 

The symbol error performance of the SCODs constructed 
in this paper (denoted as RZCOD in the plots which means 



COD with Reduced number of Zeros) for 8, 16, 32 and 64 
antennas are compared with that of well-known COD (denoted 
as SCOD) of same order in Fig.Q]and Fig.[3]under peak power 
constraint. Similarly, Fig. [2] and Fig. [4] compare the corre- 
sponding codes under average power constraint. The average 
power constraint performance of RZCOD matches with that 
of the comparable SCOD, while the RZCOD performs better 
than the corresponding SCOD under peak power constraint as 
seen in the figures. We also observe that the performance of 
our code for 8 transmit antennas matches with that of the code 
Gy (denoted as Yuen(8)) constructed by Yuen et al. 

VI. Discussion 

We have given construction for rate s ^r SCODs for 2 a 
antennas, for all values of a, with lesser number of zero entries 
than the known constructions. When a + 1 is a power of 2, 
our construction gives SCODs with no zero entries. This case 
alone generalizes the constructions in 0, (8|, |9), iflOl which 
are only for 8 antennas. Some of the possible directions for 
further research are listed below: 

• For arbitrary values of a the fraction of zero entries in our 
codes is 1 - ^tl 2 Lio92( #r )J . We conjecture that SCODs 
with smaller fraction of zero entries do not exist. It will be 
an interesting direction to pursue to settle this conjecture. 

• Several designs including CODs have been found useful 
in systems exploiting cooperative diversity. It will be 
interesting to investigate the suitability of the codes of 
this paper for cooperative diversity. 

• We have exploited the combinatorial structure of the rows 
of the design in (HJ to obtain the codes with low PAPR. 
The interrelationship between AODs and our codes is an 
important direction to pursue. 
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Fig. 1. The performance of the RZCODs and SCODs for 8, 16 and 
32 transmit antennas and the code given in Yuen et al for 8 transmit 
antennas using QAM modulation. 
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Fig. 3. The performance of the RZCODs and SCODs for 64 transmit 
antennas using QAM modulation. 
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Fig. 5. The premultiplying matrix Q( 5 ) for 32 antennas 
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